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ABSTRACT 

By employing a simple semi-analytical star formation model where the formation rates of Population (Pop) 
I/II and III stars can be calculated, respectively, we account for the number distribution of gamma-ray bursts 
(GRBs) with high pseudo-redshifts that was derived from an empirical luminosity-indictor relationship. It is 
suggested that a considerable number of Pop III GRBs could exist in the present sample of Swift GRBs. By 
further combining the implication for the star formation history from the optical depth of the CMB photons, 
it is also suggested that only a very small fraction (~ 0.6%) of Pop III GRBs could have triggered the Swift 
BAT. These results could provide an useful basis for estimating future detectability of Pop III stars and their 
produced transient phenomena. 

Subject headings: gamma-ray burst: general, reionization, first stars 


1. INTRODUCTION 

Gamma-ray bursts (GRBs), especially the major class of a 
duration longer than 2 seconds, are the most luminous objects 
in the universe, which makes them detectable even at the edge 
of the universe. On one hand, the reported highest redshift of 
GRBs so far is up to z ~ 9.4 for GRB 090429B (Cucchiara 
et al. 2011). On the other hand, some theoretical models 
predict that GRBs could be detected at much more distance 
(Band 1993; Bromm & Loeb 2006; de Souza et al. 2011). 
Therefore, it is widely suggested that long GRBs can be used 
as a cosmological tool to probe the early universe. 

Long GRBs are also known to be associated with the col¬ 
lapse of massive stars (Galama et al. 1998; Stanek et al. 2003; 
Hjorth et al. 2003), which indicates that the GRB event rates 
could trace the cosmic star formation history either unbiasedly 
(e.g.. Chary et al. 2007) or, more probably, with an additional 
evolution effect (Daigne et al. 2006; Ksitler 2008; Salvaterra 
2009; Campisi et al 2010). Thanks to the launch of the Swift 
spacecraft, the number of GRBs with a measured redshift has 
been increasing rapidly during the past decade. This makes it 
possible to clarify the connection between the GRB numbers 
and the star formation rates (SFRs; e.g., Kistler et al. 2008; 
Cao et al. 2011; Wang & Dai 2011; Tan et al. 2013) and fi¬ 
nally to infer the SFRs at high redshifts (Totani 1997; Wijers 
et al. 1998; Lamb & Reichart 2000; Porciani & Madau 2001; 
Murakami & Yonetoku 2005; Yiiksel et al. 2008; Kistler et 
al. 2009; Wang & Dai 2009; Ishida et al. 2011; Wang et 
al. 2013) where a direct SFR-measurement is extremely dif¬ 
ficult. Nevertheless, such attempts could somewhat be dis¬ 
turbed/hindered by the limited number of measured redshifts 
and, in particular, by some inevitable observational selection 
effects (Guetta & Piran 2007; Cao et al. 2011). 

As an alternative and complementary method, in Tan et 
al. (2013), we proposed to estimate pseudo-redshifts for the 
Swift GRBs with an empirical relationship between the spec¬ 
tral peak energy and the peak luminosity. Consequently, 498 
pseudo-redshifts up to z ~ 30 were derived from the Butler’s 
GRB catalog where the spectral peak energies are roughly 


estimated with Bayesian statistics rather than observed. Al¬ 
though the L - E p relationship is not so tight and the obtained 
pseudo-redshifts are not so reliable individually, a statisti¬ 
cal study on the pseudo-luminosity distributions of the GRBs 
can in principle be implemented, especially worthy of men¬ 
tion, for different narrow redshift ranges. As a result, the 
luminosity function of the GRBs was found to be strongly 
redshift-dependent, for redshifts z < 3.5 and a previously- 
determined star formation history. From the pseudo-redshift 
sample of Tan et al. (2013), we can find that there are 38 
GRBs whose pseudo-redshifts are higher than 10, where 23 
GRBs are above the lower luminosity cut off of Swift satellite. 
This statistically indicates that a remarkable number of high- 
redshift GRBs may exist in the present Swift GRB sample, 
most of which, however, evaded direct redshift measurements 
due to redshift selection effects. 

As widely considered, for redshifts z > 10, high mass Popu¬ 
lation (Pop) III stars are probably more dominant than normal 
Pop I/II stars (Bromm et al. 2002; Abel et al. 2002). Pop 
III stars were born in metal free gas and their deaths started 
the metal enrichment of the intergalactic medium (IGM) via 
supernova feedback, which subsequently lead to the forma¬ 
tion of Pop I/II stars (Ostriker & Gnedin 1996; Madau et 
al. 2001; Bromm et al. 2003; Furlanetto & Loeb 2003). 
The high mass and zero-metallicity of Pop III stars make 
them easily to produce a collapsar GRB (Hirschi 2007), 
whose isotropically-equivalent energy could be as high as 
~ 10 55 ~ 57 ergs (Meszaros & Rees 2010; Suwa & Ioka 2011). 
Additionally, as the first generation stars in the universe. Pop 
III stars also turn on the cosmic reionization by emitting ultra¬ 
violet photons, and so the reionization history should strongly 
depend on the formation history of Pop III stars (Furlanetto & 
Loeb 2005; Barkana 2006; Robertson et al. 2010). 

Therefore, in this paper we try to use the GRBs with high 
pseudo-redshifts (z > 3.5) to probe the SFRs up to the begin¬ 
ning of the reionization era as well as the detection efficiency 
of Pop III GRBs, by combining the constraint from the optical 
depth of the cosmic microwave background (CMB) photons. 
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TABLE 1 

The GRB numbers in different redshift ranges. 


Redshift 


N* 

p 

0-0.5 

10 

28 

0.5 - 1 

30 

74 

1 - 1.5 

22 

96 

1.5-2 

23 

65 

2 - 2.5 

20 

37 

2.5 - 3 

19 

37 

3-3.5 

14 

24 

3.5 -4.5 

12 

24 

4.5 - 5.5 

5 

14 

5.5 - 6.5 

2 

7 

6.5 - 8 

2 

9 

8-10 

1 

11 

10 - 14 

0 

12 

14 - 18 

0 

5 

18 - 22 

0 

2 

22 - 26 

0 

1 

26 - 30 

0 
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'GRBs of observed redshift. 
T’jRBs of a pseudo-redshift. 


In the next section, we will estimate the high-redshift SFRs 
by a usual method where the contribution to GRBs by Pop 
III stars is ignored, which could lead to an unacceptable re¬ 
sult. In contrast, in section 3 we attempt to understand the 
GRB numbers by considering the GRB productions of both 
Pop I/II and III stars, where a simple semi-analytical model 
for the collapse of dark matter halos is employed. Finally, the 
results are given and discussed in section 4. 


2. SFR-DETERMINATION WITHOUT POP III STARS 

Table 1 lists the numbers of GRBs with observational or 
pseudo-redshifts for different redshift ranges, which are taken 
from Tan et al. (2013). Here only GRBs with luminosities 
higher than the lower luminosity cut off of Swift are taken into 
account. The difference between these two sets of numbers 
arises from the selection effects of redshift measurements that 
are, however, difficult to be described theoretically. In this 
paper, we will only pay attention to the numbers of pseudo- 
redshifts in order to probe the high-redshift ranges and avoid 
the redshift selection effects. 

In view of the model, for a given GRB luminosity function 
ffi-(L) and GRB event rate R(z), the expected number of GRBs 
within redshift range z,\ < z < Zi can be calculated by 


^kl,Z2] — 


—r r r 

4?r J z j Jl Ic (z) 


O z (L)R(z)dL 


dV{z) 

1 +z ’ 


( 1 ) 


where AQ/47 r ~ 0.1 is the field view of the Burst Alert Tele¬ 
scope (BAT) on board Swift, T ~ 7yr is the observational pe¬ 
riod, L] c (z) is an adopted lower cutoff of the GRB luminosity 
corresponding to the selected data [see Equation (3) in Tan et 
al. (2013) and explanation therein], and dV(z) is the comov¬ 
ing cosmic volume. 

From the GRBs having a pseudo-redshift z < 3.5 and the 
corresponding star formation history 


(1 + z) 3 ' 44 , z < 0.97, 

(1 + z)°, 0.97 <z< 3.5, 



Fig. 1.— SFRs inferred from the GRBs with pseudo-redshifts > 3.5 (solid 
circles), where the contribution to the GRBs by Pop III stars is ignored. The 
crosses display the SFRs according to Hopkins & Beacom (2006) and the 
open diamonds correspond to the GRB-inferred SFRs given by Kistler et al. 
(2009). The solid (or dashed) line represents a power-law star formation 
history, which is required by the CMB optical depth for N y = 4000 and 

/esc = 10% (or /esc = 2%). 


Tan et al. (2013) derived a GRB luminosity function by: 


0 ; (L) = 


1 

4.5L b (z) 



L < L b (z), 
L > L b (z), 


(3) 


where the break luminosity reads L b (z) = 1.2 x 10 51 (1 + 
z) 2 erg s' 1 , and a GRB event rate as 

R(z) = f B p(z)p*(z), (4) 

where the proportionality coefficient f B p(z) = 2.4 x 10~ 8 (1 + 
zY'M q x with f B being the beaming factor of the GRB jets and 
p(z) representing the GRB-production efficiency of the stars. 

As the most straightforward assumption in this paper, we 
assume that Equations (3) and (4) are valid for all redshifts, 
which could be correct if all GRB progenitors belong to the 
same type of stars (i.e.. Pop I/II stars). Then following Kistler 
et al. (2008, 2009), the SFRs at high redshifts can be esti¬ 
mated by the following equation 


P*(z ) 


%. 3 .aGf L^WpizidL^ 


(5) 


Substituting the numbers of GRBs with pseudo-redshifts for 
different redshift ranges into the above equation, we derive the 
high-redshift (z > 3.5) SFRs as shown by the solid circles in 
Figure 1, where a big bump appears surprisingly in the cosmic 
star formation history within the redshift range 8 < z < 30. 

By considering the connection between the star formation 
and cosmic reionization, the SFRs can also be constrained 
by the CMB optical depth r = 0.088 + 0.015 measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP) exper¬ 
iment. Supposing a simple power-law high-redshift star for¬ 
mation history asp*(z) oc (1 +z) _ff forz > 3.5, Yu et al. (2012) 
derived a constraint on the index a as 


cr = 2.18 lgWy/e SC ) - 3.89, (6) 

where N y is the number of ionizing ultraviolet photons re¬ 
leased per baryon of the stars and / esc is the fraction of these 
photons escaping from the stars. For a typical value of N y ~ 
4000 for a Salpeter stellar initial mass function and a metallic - 
ity 0.05Z o (e.g. Barkana 2001), we can get a — 1.78 and 0.26, 
corresponding to a typical escaping fraction / esc ~ 10% and 


p*(z) or 
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a very small one / esc ~ 2%, respectively. Such requirements 
on the SFRs are shown by the lines in Figure 1, which are, 
however, significantly lower than the GRB-inferred SFRs. In 
other words, the SFRs simply inferred from the GRB num¬ 
bers would lead to an overionized universe, unless the pho¬ 
ton escaping fraction is unacceptably smalQ This indicates 
that something wrong must appear in the above processes of 
determining the SFRs by GRBs. The most probable reason 
could be that it is misleading to assume Equations (3) and (4) 
to be valid for all GRBs. In fact, for redshifts >10, Pop III 
stars could play a dominant role in producing GRBs. Due to 
the unique properties of Pop III stars, the formation history 
of Pop III GRBs is probably very different from that of Pop 
I/II GRBs. In such a case, one may worry whether the Pop 
III GRBs can obey the same L - E p relationship as Pop I/II 
GRBs. In our opinion, the L - E p relationship is probably de¬ 
termined by the emission mechanism of the GRB jets, which 
could not be sensitive to the progenitor properties. 

3. CONFRONTING THE GRB NUMBERS WITH A STAR 
FORMATION MODEL 

3.1. Semi-analytical star formation model 

In the hierarchical formation model, star formation takes 
place during the collapsing and merging of dark matter ha¬ 
los. A straightforward semi-analytical approach for the abun¬ 
dance of dark halos was first given by the Press-Schechter 
(PS) formalism (Press & Schechter 1974), which was subse¬ 
quently improved by Sheth & Tormen (1999). By using the 
Sheth-Tormen mass function (ST) nsr{M, z) of dark halos, the 
collapse fraction of mass available for star formation can be 
calculated by 

/cot(z)= I n ST (M, z)—dM, (7) 

JM m i„ P 

where p is the mean density of the universe and M m j n is the 
minimum mass of the halos below which the halos can not 
collapse. By requiring the virial temperature 7\. lr of the halos 
to be higher than the temperature permitted by possible cool¬ 
ing channels, the minimum halo mass can be determined to 

M mm ~ 10 8 M o (r vir /10 4 K) 3/2 [(1 +z)/10]- 3/2 . Here we typi¬ 
cally assume the minimum halo mass corresponds to the virial 
temperature of 7\. lr = 10 4 both for Pop I/II and Pop III stars, 
above which atomic hydrogen cooling is effective. 

The star forming halos (“galaxies”) could launch a wind of 
metal-enriched gas with an initial speed of vo by producing a 
large number of supernovae. The well-known Sedov (1959) 
solution predicts that the galactic wind could travel a comov¬ 
ing distance of R w = v w t = (3E w /2nfbp) 1/5 t 2/5 , where v w is 
the average wind speed while traveling in the apace, p is the 
density, E„ = is the total energy of the galac¬ 

tic wind with f, = 0.167 being the mass fraction of baryonic 
matter, f, the star formation efficiency, and / w the energy frac¬ 
tion that goes into the wind. Here we have assumed that each 
wind has propagated for half of the age of the universe with 
t ~ l/377(z). Therefore, the volume-filling fraction of the 
metal-enriched gas of a collapsing halo can be expressed by 
rj(M,z) = (Rvi/Rh) 3 , where R/ t = (3M/7207rp)^ 3 is the radius 

1 For an escaping fraction less than 2%, the CMB optical depth would 
require much higher SFRs (a < 0.26). However, in such cases, the cosmic 
reionization can only be complete at redshifts < 6, which is inconsistent with 
the result from the probes to the Gunn-Peterson trough (Ly-o' absorption) 
toward high-redshift quasars and galaxies (Yu et al. 2012). 



1+z 

Fig. 2.— Fitting to the number distribution of the GRBs with pseudo- 
redshifts (grey histograms) by combining the contributions to GRB produc¬ 
tions from Pop I/II (dash-dotted line) and III (solid line) stars. 


of the whole halo (e.g. Johnson 2010). Then, by consider¬ 
ing the mass distribution of the halos, the total fraction of the 
cosmic volume that is metal-enriched by the galactic winds 
can be written as (Furlanetto & Loeb 2005; Greif & Bromm 
2006) 

Qeiz) = [1 + fth(z)] f r](M,z)nsj(M,z)—dM, (8) 
Jm m „ P 

where fhhiz) is the linear correlation function between two 
halos with mass M at a comoving distance R: 

fth (M, R, z) = b 2 (M , z)f mm (R), (9) 


here b(M , z)is the linear bias of a halo and f mm (K) is the mass 
correlation function (e.g. Mo & White 2002) with R being the 
average wind size in the comoving unites. The average bias 
of the enriched regions is given by. 


C dMMrj(M, z)b(M, z)«st 

Hz) =-—-- 

I dMMnsr 


GO) 


By assuming that the host galaxies were randomly distributed, 
we could get the enrichment probability p e = 1 -exp [~Q e (z)]. 

Since the halos hosting zero-metallicity Pop III stars should 
not be located in the wind radius of old galaxies, the SFR of 
Pop III stars at redshift z can be roughly estimated by 


P*\z) = fl"fbP( ] ~Pe ) 


d.f a , dz 
dz dt ’ 


( 11 ) 


where /J 11 is the fraction of baryonic mass that goes into Pop 
III stars. Simultaneously, the SFR of Pop I/II stars read 

p'i\z)=fr/ h p P ^ d i, G2) 

dz dt 

where / 4/H is the baryonic fraction for Pop I/II star formation. 


3.2. Fitting to the GRB number distribution 

Different from the treatment in Section 2, here we assume 
the GRB luminosity function presented in Eq. (3) and the 
connection between the GRB event rate and SFR in Eq. (4) 
are only valid for the GRBs produced by Pop I/II stars. Instead 
of Eq. (1), we calculate the number of Pop III GRBs within 
redshift range j| < z < z.z by the following formula: 


N 


in 

[Zl,Z2] 


AQ r Zl -TIT dV(z) 

= — TfJ R m (z)-r—dz 


1 +z 


(13) 
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with R m (z) = ftip m p ] * l (z), where the parameter f repre¬ 
sents the detection efficiency of Pop III GRBs (see discus¬ 
sions at the end of Section 3.3) and the GRB-production ef¬ 
ficiency of Pop III stars p m is considered to mainly arise 
from the mass requirement on the progenitors. The masses 
of Pop III stars are usually considered to be higher than 
100M o (Bromm et al. 2002; Abel et al. 2002), but a much 
lower mass as ~ 10M o was also suggested by some re¬ 
cent more elaborate researches where the feedback effects 
are taken into account simultaneously (Greif & Bromm 2006; 
de Souza 2011; Hosokawa et al. 2011, 2012; Stacy et al. 
2012). In such a case, there could be only a small fraction 
of Pop III stars available for GRB production, because the 
GRB production could require the progenitor’s mass higher 
than ~ 40M o (Belczynski et al. 2007; Tornatore et al. 
2007; Marassi 2009). In addition, it is widely suggested that 
Pop III stars with masses 140M o < m < 260M o would 
be completely disrupted by a supernova explosion due to 
the pair-creation instability (PISN; Heger & Woosley 2002), 
so they also can not contribute to GRBs. Finally, by in¬ 
voking the Salpeter mass function f{m), the GRB produc¬ 
tion efficiency of Pop III stars can be estimated by p m = 


140 M 0 
'40 M 0 


X 


<p(m)dm + /“ 4>(m)dm 


u 


10 3 M q 

10M o 


M(f>(m)dm — 4 x 


ur 3 M-'. 

Substituting Eq. (12) into (1) and Eq. (11) into (13), we can 
calculate the numbers of Pop I/II and III GRBs, respectively, 
for specific values of model parameters. By confronting these 
model-predicted GRB numbers with the numbers listed in Ta¬ 
ble 1, we can in principle obtain an observational constraint on 
the model parameters. We firstly take typical values for some 
parameters according to some previous works as 7\. lr = 10 4 
K for atomic hydrogen cooling (Haiman et al. 2000; Furlan- 
etto & Loeb 2005), f w = 0.1 (e.g., Scannapieco & Broad- 
hursts 2001; Springel & Hernquist 2003;), and v'o = 200km/s 
(Fabian et al. 1980; Eymeren et al. 2007). For the remain¬ 
ing more uncertain parameters fl /n , /J n , and f, their best-fit 
values are constrained to f! /n =0.16 and f} u f - 6 x 10 -6 by 
fitting to the GRB number distribution without error bars. The 
corresponding fitting is shown in Figure 2 by the lines. The 
value of f* 11 could be well consistent with some previous an¬ 
alytical estimations and simulations (Greif & Bromm 2006; 
Marassi et al. 2009). Qualitatively, our result indicates that 

(1) the formation history of Pop I/II stars with typical values 
of parameters is basically favored by the GRB numbers, and 

(2) a considerable number of Pop III GRBs could exist in the 
present Swift GRBs although the degeneracy between f and 
ff' has not been removed. 


3.3. Cosmic reionization 

In order to remove the degeneracy between the parameters 
f and /„ m , we further invoke an implication for the star for¬ 
mation history from the CMB optical depth. The evolution 
of the cosmic reionization denoted by x = «h ii/«h can be 
determined by the following equation 


dx 

dz 


-— y —— - a B C(z)( 1 + z) 3 (l + y)u„x 
(i + y) H H 


dr 

dz 


(14) 


with a rate of ionizing ultraviolet photons escaping from stars 
into IGM hy(z) = Nyf esc p»(z)/ni B , where m B is the mass of 
a baryon. y = 0.08 by assuming that the helium was only 
once ionized, - 1.9 x 10~ 7 cm 3 is the local number den- 



Fig. 3.— High-redshift star formation history of Pop I/II (dash-dotted line) 
and III stars (solid line), which are constrained by the GRB numbers and the 
CMB optical depth together. For comparisons, some simulation results of the 
formation history of Pop III stars are also presented by dashed lines, such as 
Tornatore et al. (2007; T07), Trenti & Stiavelli (2009; TS09), and Kulkami 
et al. (2013; K13). The dotted lines with arrows give an upper limit on the 
SFRs of Pop III stars, which is taken from Inoue et al. (2013). 

sity of hydrogen, a B = 2.6 x 1CF 13 cm 3 s' 1 is the recombina¬ 
tion coefficient for electron temperature of about 1 0 4 K, and 
C(z) = 2.9 [(1 + z)/6] 11 (or = 2.9) for z > 5 (or z < 5) is 
the clumping factor of ionized gas. With a given reionization 
history x(z), the CMB optical depth can be calculated by inte¬ 
grating the electron density times the Thomson cross section 
along proper length as 

t =-(l + y)o- T n° H c f (1 + zfx(z)^-dz. (15) 

Jo dz 

Here we define the upper limit of the integral to be a moderate 
value with Zh ~ 30, because the CMB optical depth could be 
mainly contributed by the electrons at relatively low redshifts 
Z <sc Zh, and also as the WMAP experiment is insensitive to 
too high redshifts (Larson et al. 2011). 

In the star formation model with the parameters constrained 
above, the remaining free parameter /J n can be determined 
to be /JH = 10 3 from the CMB optical depth r = 0.088, 
where we take fUf = 6%, (V 3/u = 4000, f™ = 30%, and 
A 311 = 3 x 10 4 (Greif & Bromm 2006). Such a value is just lo¬ 
cated within its theoretically-expected range of ~ 10 6 - 10 3 
(Greif & Bromm 2006; Marassi et al. 2009). At the same 
time, the detection efficiency of Pop III GRBs is revealed to 
be f = 0.6%. Such a low efficiency could be caused by var¬ 
ious reasons. For example, the luminosity selection with an 
unknown luminosity function of Pop III GRBs should be in¬ 
cluded in the parameter J, and some Pop III GRBs might oc¬ 
cur in extremely dense accretion envelopes that suppress the 
GRB luminosity at early times. Moreover, the high gamma- 
ray variability of the GRBs could be significantly smoothed 
by the high redshifts, which could lead to facilities such as 
Swift BAT not to be triggered by the gamma-ray emission. 

4. RESULTS AND DISCUSSIONS 

The star formation histories of Pop I/II and III stars with the 
obtained parameter values are presented in Figure 3. Firstly, 
the SFRs at redshift range 5 < z < 10 are suggested to be 
slightly higher than the ones given by Kistler et al. (2009), 
because here an evolving GRB luminosity function is adopted 
for Pop I/II GRBs and a larger GRB number with pseudo red- 
shift is used here. Secondly, the formation history of Pop III 
stars is demonstrated to be basically consistent with some pre¬ 
vious simulations and an observational upper limit is derived 
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from the gamma-ray attenuations. To be specific, the peak 
of the SFRs of Pop III stars appears from our analysis to be 
at z ~ 11 and the formation of Pop III stars could continue 
down to as low as z ~ 5. The transition from the formation 
of Pop III stars to Pop I/II stars takes place at z ~ 18. Fi¬ 
nally, although the simple star formation model can account 
for the GRB numbers for a wide redshift range as z < 25, the 
GRB number for 26 < z < 30 still significantly exceeds the 
model prediction. Such an obvious excess may indicate some 
other GRB origins in the early universe, e.g., superconducting 
cosmic strings (Cheng et al. 2010). In this paper, we do not 
provide a precise regression analysis of the data, by consider¬ 
ing that the large intrinsic uncertainty of the present pseudo- 
redshift data makes it unnecessary to find a high-precision fit¬ 
ting. In the future, a more accurate method to get the GRB 
pseudo-redshift is expected. 

The determination of the SFRs of Pop III stars is of fun¬ 
damental importance for estimating the detection efficiencies 
of these stars and their produced transient phenomena with 
some future facilities. For example, the upcoming James 
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